Micromagnetic theory of spin relaxation and ferromagnetic resonance in 
multilayered metallic films 

Sergey Bastrukov, 1a) Jun Yong Khoo, 1 Boris Lukiyanchuk, 1 and Irina Molodtsova 2 

1) Data Storage Institute (DSI), Agency for Science, Technology and Research (A*STAR), 
5 Engineering Drive 1, 1 17608 Singapore 

2) Laboratory of Informational Technologies, Joint Institute for Nuclear Research, 
141980 Dubna, Russia 

(Dated: 22 October 2012) 

Spin relaxation in the ultrathin metallic films of stacked microelectronic devices is inves- 
tigated on the basis of a modified Landau-Lifshitz equation of micromagnetic dynamics in 
which the damping torque is treated as originating from the coupling between precessing 
magnetization-vector and the introduced stress-tensors of intrinsic and extrinsic magnetic 
anisotropy. Particular attention is given to the time of exponential relaxation and ferro- 
magnetic resonance linewidth which are derived in analytic form from the equation of 
magnetization energy loss and Gabor uncertainty relation between the full-width-at-half- 
maximum in resonance-shaped line and lifetime of resonance excitation. The potential of 
developed theory is briefly discussed in the context of recent measurements. 
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I. INTRODUCTION 



An understanding of relaxation processes in ultrathin films of ferromagnetic metals is crucial 
to the design and construction of microelectronic devices 1,2 , like magnetic random access memory 
(MRAM) and spatial light modulators (SLM). The main source of information about relaxation 
processes has been and still is the ferromagnetic resonance (FMR) measurements aimed at reveal- 
ing the frequency dependence of the full-width-at-half-maximum in FMR spectral line 3 ' 4 . Tradi- 
tionally, the results of these experiments are treated within the framework of the phenomenological 
Landau-Lifshitz-Gilbert model", describing the FMR response in terms of uniform precession of 
magnetization M(i) with the preserved in time magnitude |M(t)| = M s where M s is the satu- 
ration magnetization. The dynamical equation governing precessional motions of M(t) can be 
conveniently written in the following general form 

M(t) = - 7/ i T(t)-R(t), (1) 

where 7 is the electronic gyromagnetic ratio and fi is the magnetic permeability of free space; 
SI units are used throughout this paper. The vector-function, T(t) = [M(t) x H], represents the 
magnetic torque that drives the free Larmor precession of M(t) about the axis of the dc magnetic 
field, H = constant, in the process of which the Zeeman magnetization energy, W m (t) = — /i H • 
M(t), is conserved: W m (t) = 0. Central to understanding the relaxation process is the relaxation 
function R(£) defining the rate of magnetization energy loss 

W m (t) = -// H • M(t) = /i H ■ R(t). (2) 

In this work we focus on Landau-Lifshitz (LL) form of this function Rj n (t) = \ in [M.(t) x [M(t) x 
H]] which provides geometrically transparent insight into the magnetization- vector motion in the 
process of aligning M with H. The material-dependent parameter X in can be thought of as de- 
scribing the strength effect of the intrinsic anisotropy on the relaxation dynamics of magnetization 
precession which is constrained by the conditions M(t) ■ M(£) = and M(£) • R(£) = 0. In this 
paper we consider an alternative micromagnetic treatment of R(t) according to which the origin 
of the damping torque responsible for spin relaxation in multilayered metallic films is attributed 
to the coupling between the uniformly precessing magnetization-vector and the stress-tensors of 
intrinsic and extrinsic magnetic anisotropy. 
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H. STRESS-TENSOR REPRESENTATION OF MICROMAGNETIC DAMPING 
TORQUE 



The equilibrium magnetic anisotropy exhibited in the easy and hard axes of magnetization di- 
rection 8 is a hallmark of ultrathin films of ferromagnetic metals. Viewing this property from the 
perspective of the macroscopic electrodynamics of magnetic continuous media 910 , it seems quite 
natural to invoke the stress-tensor description of magnetic anisotropy, namely, in terms of sym- 
metric tensors of magnetic-field-dependent stresses. In so doing we adopt the following definition 
of stress-tensor of intrinsic anisotropy (generic to both monolithic and multilayered ferromag- 
netic films) and stress-tensor of extrinsic anisotropy (arising from impurities and imperfections 
of the film crystalline lattice) 

< = \ l(M n H n )5 kl - (M t H k + M k H t )} , (3) 
<r% = | [H 2 5 kl - (HiH k + H k Hi)] , (4) 

where 5 ik is the Kronecker symbol and fi stands for the effective magnetic permeability which is 
derived from the magnetization curve according to the rule 11 : /x = AB/ AH. Hereafter H refers 
to the total (applied and internal effective) field. It is worth noting that the above stress-tensor 
description of intrinsic and extrinsic magnetic anisotropy is consistent with the definition of the 
energy density of magnetic field stored in a ferromagnetic film 12 

m=^B-H, B = ii(H + M) (5) 

in the sense that the relation between the stress-tensor of combined intrinsic and extrinsic magnetic 
anisotropy, a ik = oj£ + afg, and the energy density u is described by 

u = Tr[a lk ] = a u = ^(MH + H 2 ), (6) 

where Tr[ai k ] stands for the trace of tensor a\ k . In what follows we focus on the effect of above 
magnetic stresses on the precessing magnetization vector whose mathematical treatment is sub- 
stantially relied on the symmetric tensor 

lik = [M 2 J lk - MiM k ], llk = 7fci , (7) 

having, in appearance, some features in common with that for isotropic magneto- striction stresses 13 
It can be verified by direct calculation that the stress-tensor representation of the intrinsic relax- 

3 



ation function is identical to the LL relaxation function 

2A 

R T = -i^Hk Mi a% = X in [Mi{M k H k ) - HtM?], (8) 
R m = \ in [M(t) x [M(t) xH]]. 

In choosing the above form of the tensor j ik we were guided by previous investigations 14 of the 
damping terms in ferro-nematic liquid crystals dealing with the tensor constructions of a similar 
form. For the extrinsic relaxation function, owing its origin to the coupling of Mi with af%, we use 
the following stress-tensor representation 

Rf = -^HkMtf* (9) 

= -\J-^§^[M t {M k H k ) - HiM?]. 

The minus sign means that extrinsic damping torque counteracts the damping torque originating 
from the intrinsic stresses. The vector form of extrinsic relaxation function (9) reads 

(M(t) • H) . , , 

Rex(t) = -Aea,- >- [M(t) X [M(t) X H]]. (10) 

As is shown below, the parameter-free frequency of the transient magnetization configuration 

to 

^ = i 7 T^TT^' u m = 7/^o^"s, w = 27r/, (11) 

provides correct physical dimension of the extrinsic damping torque and proper account for the 
empirical dependence of the FMR linewidth AH upon the resonance frequency /. Making use of 
argument of physical dimension it is easy to show that the material-dependent parameters X in > 
and X ex > (measuring strength of intrinsic and extrinsic stresses on the relaxation process in 
multilayered film) can be represented in terms of dimensionless damping constants a and (3 (whose 
magnitudes are deduced from the empirical frequency dependence of FMR linewidth) as follows 

The net outcome of the above outlined procedure of computing the combined (intrinsic plus ex- 
trinsic) damping torque 

R- = R-in + Rex (13) 

(M(t) • H) 



Xin X ex 



[M(t) x [M(i) x H]], 
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entering the basic equation of micromagnetic dynamics (1) is the following Modified Landau- 
Lifshitz (MLL) equation 



M 



7yU [M x H 
7/^o 



(14) 



a 



7/^o 
Mo 



(M-H) 



[M x [M x H]], 



which obeys all constraints of the canonical LL equation. One sees that unlike the linear-in- 
magnetic-field intrinsic damping torque, the extrinsic damping torque is described by quadratic- 
in-magnetic-field relaxation function. At this point it seems noteworthy that the need in allowing 
for the quadratic-in-iJ damping terms has been discussed long ago 15 . The above scheme can be 
regarded, therefore, as a development of this line of theoretical investigations. In terms of the unit 
vector of magnetization m(t) = M(t)/M s and Larmor frequency uj = 7/ioH the last equation 
can be converted to (see 16 for comparison) 



m 



\uj x m 



a — (3 



m ■ us] 



n 



[m x [m x oj]}. 



(15) 



It can be seen that the obtained MLL equations (14) and (15) are reduced to the standard LL 
equation when the effect of extrinsic stresses is ignored (i.e. (3 = 0). 



III. VARIATION METHOD OF COMPUTING RELAXATION TIME AND FMR 
LINEWIDTH 

The relaxation time is amongst the primary targets of current FMR experiments. In this section, 
we present variational method of analytic computation of the FMR linewidth which is quite dif- 
ferent from the well-known solution of the susceptibility solution of LL equation 7 . At the base of 
the variation method under consideration lies the equation of the magnetization energy loss from 
which the exponential relaxation time r as a function of FMR frequency / = u;/(27r) is derived. 
The FMR linewidth, Aco = 7/ioAif , is computed from the well-known Gabor uncertainty relation 
(e.g. 17 ,Sec.ll.2) 

Awr=l (16) 

between the full-width-at-half-maximum in the resonance- shaped spectral line Aco and lifetime r 
of resonance excitation. 
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A. FMR linewidth caused by intrinsic damping torque 



For the former we consider relaxation process brought about by intrinsic damping torque. Our 
approach is based on the observation that the equation of magnetization energy loss in the process 
of a uniform precession of magnetization in a dc magnetic field 

™=. = -, HM,^M = ^mt), (17) 

K = -^7*Mi°£, (18) 

is reduced to the equation for the cosine function u(t) = cos#(£) of angle 9(t) between M(£) and 
H, namely 

^ = -auj[u 2 (t) - 1], u = 1 ^H. (19) 

The right hand side of (19) suggests that there are two equilibrium configurations, namely, with 
u(0) = 1 corresponding to M || H and u(ir) = —I corresponding to M t! H. The stability 
of these configurations can be assessed by the standard procedure of introducing small- amplitude 
deviations 5u[t) from the equilibrium values u(0) = u = ±1. On substituting u(t) = u + Su(i), 
into (19) with u — 1 and retaining first order terms in Su(t) we obtain equations describing 
exponential relaxation of magnetization to the state of stable magnetic equilibrium: 

= -(2au)5u -> Suit) = 5 M (0)e-' /r , (20) 

dt 

r _1 = 2atu, u} = 7/io-ff- (21) 

The second stationary state, with u = —1, is unstable, since in this case the resultant linearized 
equation, 5u = (2au)5u, having the solution, Su(t) = 5w(0)e*/ T , describes a non-physical behav- 
ior of Su as the time is increased. Inserting (21) in (16), we arrive at the basic prediction of the 
standard micromagnetic model 

AH(f) = Af, A = —. (22) 

This last equation provides a basis for discussion of empirical linewidth-frequency dependence 
A.H exp = AH + AH exp (f) with AH exp (f) = A exp f. Central to such a discussion is the identi- 
fication of theoretical and experimental linewidths, AH(f) = AH exp (f), from which the magni- 
tude of a is deduced and applied to (21) for obtaining numerical estimates of the relaxation time 
r. 
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B. FMR linewidth caused by both intrinsic and extrinsic damping torques 



In this case the starting point is the equation of magnetization energy loss with the combined 
relaxation function 



dWr, 



dt 



—ji HM t 



Ri 



d(cos0(t)) 
It 



VoHiRi(t), 



2A 



(23) 
(24) 



which after some algebra is converted into equation for u having the form 



du{t) 
dt 



aid — f3—u(t) 



(u 2 (t) 



(25) 



The right hand side of this equation suggests that there are three stationary state characterized by 



u (6 = 0) = ±1, u (6 = e M ) 



0u' 



(26) 



Applying to (25) the standard linearization procedure u{t) = u + Su(t) in (25) one finds that 
resultant equation is equivalent to the equation of exponential relaxation, 5ii = — r" 1 ^, if and 
only if the parameter 



T 



UJ 2 \ UJ 2 



(27) 



is a positive constant. It is easy to see that this is the case for uq(9 = 0) = 1 and uq{6 = 9m) 
given by rightmost of equations (26). This latter u corresponds to a quasi- stationary transient 
configuration of precessing magnetization owing its existence to the coupling of magnetization 
with extrinsic stresses of magnetic anisotropy. The state with uq — — 1, is unstable. For the total 
relaxation time r _1 = T^iO = 0) + r _1 (6 ) = M ) and the FMR linewidth (following from Gabor 
uncertainty relation AH = [7/i r] -1 ) we obtain 

,2 



2aw-/3^-(l + cos 2 9 M ) 



AH 



2cu 

4vr/ 
7/^o 



a 



a 




(28) 
(29) 



'l + cos 2 6 



Mj 



It is worth emphasizing that the expounded micromagnetic mechanism of the magnetization 
precession damping (due to magnetization-stress coupling) presumes that the process of spin- 
relaxation is not accompanied by generation of spin- waves (magnons), because the magnetization 
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FIG. 1. The FMR linewidth AH and relaxation time r as functions of the FMR frequency /, computed on 
the basis of the standard and modified in the present work Landau-Lifshitz equation. 

M is regarded as a spatially-uniform vector across the multilayered film. At this point the con- 
sidered regime of the magnon-free spin relaxation (in which the wave vector of spin wave k = 0) 
is quite different from spin relaxation caused by two-magnon scattering 18 . The most conspicuous 
feature of this (substantially macroscopic) mechanism, responsible for the non-linear frequency 
dependence of FMR linewidth, is the transient magnetization configuration owing its existence to 
the extrinsic stresses generic to the multilayered films. Such a configuration is absent in perfect 
monolayered films (without impurities and defects of crystalline lattice) of pure ferromagnetic 
metals (Ni, Co, Fe) whose ferromagnetic properties are dominated by intrinsic stresses of mag- 
netic anisotropy. 

IV. DISCUSSION AND SUMMARY 

In approaching the interpretation of FMR measurements in terms of presented theory, in the 
remainder of this work, we focus on a case of in-plane configuration (M ft H) which is of 
particular interest in connection with the recent discovery of non-linear frequency dependence of 
FMR linewidth 419 20 . In this case, the last equation for the FMR linewidth takes the form 

AH = [a — P(l - ( 7 /x M s /27r/) 1 / 2 )] . (30) 

To illuminate the difference between predictions of the standard and modified LL models, in 
Fig.l we plot t and AH as functions of the FMR frequency / computed with the pointed out 
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FIG. 2. Theoretical fit, equation (30), of the empirical non-linear frequency dependence of FMR linewidth 
detected in the FMR measurements 19 on multilayered metallic nanostructures Pd/Fe/GaAs. 
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FIG 3. Same as Fig.2, but for the thin film of Fe/V measured in 



parameters of a and f3. In computation based on the standard micromagnetic model, equation 
(22), we have used one and the same value of parameter a as in 3 reporting the FMR measurements 
on ultrathin films of Permalloy. The presented in Fig. 1 values of a and (3 have been deduced from 
fitting, equation (30), of the non-linear frequency dependence of FMR linewidth discovered in the 
FMR measurements 19 . The result of this fit is shown in Fig.2. In Fig. 3, we plot our fit of the 
FMR linewidth measurements 20 on multilayered samples of Fe/V. A more detailed discussion of 
consequences of considered micromagnetic mechanism of spin relaxation will be the subject of 
forthcoming article. 
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